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ABSTRACT 


A technique for calculating the near fields of reflector 
antennas using the geometrical theory of diffraction is 
developed* Computed patterns for two antennas are compared 
with patterns obtained by other techniques and experimental 
data*. The GTD is proved to be a useful tool in obtaining 
near field patterns of antennas* 
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I, Introduction 


The study of the near fields of antennas has several 
practical applications. Siting high power radars demands a 
knowledge of the near field pattern of the system so that 
hazards to personnel and flammable or explosive material may 
be calculated and satisfactory safety margins established. 

Since a radome is, almost exclusively, located in the near 
field of the antenna it houses, its effect on the pointing 
accuracy and beam shape of the antenna depends on the near 
field distribution at the surface of the radome. Thus a de- 
sign minimizing pointing errors and pattern degradation must 
be based on a knowledge of the near fields of the antenna. 

With the advent of millimeter wave experiments (16 GHz - 90 GHz) 
the near field. of an antenna may extend a large distance from 
the antenna. For example, the Fresnel zone of a thirty-foot 

diameter antenna operating at 32 GHz extends to 58,000 feet 

2d2 

as obtained from the far field condition atmos- 

pheric disturbances affect ground to satellite communication 
entirely through interaction with the near field of the 
antenna. The simplifying approximations to the radiation 
integrals which enable simple calculation of the far field 
patterns of large antennas are not applicable in the near field. 
Solution of the exact radiation integrals in the near field re- 
quires a sizeable investment of computation time and machinery 
for all but the smallest antennas* The computation time and 
the size of the machine required increases sharply with the 
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size of the antenna. The geometrical theory of diffraction 
provides an attractive alternative to the classical theory for 
the solution of near field problems. The techniques used in 
this study reduce the area integration required in the exact 
method to, at worst, a line integration and, at best a sum- 
mation of a few discrete terms. 

The problem is restricted here to the determination of 
the forvsrard near fields of reflector antennae in a cylindrical 
region with diameter twice the antenna diameter using the 
geometrical theory of diffraction. This region is treated 
since it represents the portion of the task that is not clearly 
included in previous solutions. The geometrical theory of 
diffraction is a logical extension of the theory of geometrical 
optics; it introduces a class of rays generated by the inter- 
action of the incident field with a discontinuity of the bound- 
ary (in our case a sharp edge.) An intuitive analysis of the 
effects of the edge shows that the diffraction products act 
to remove the discontinuity in the fields at the boundary be- 
tween the illuminated, and shadowed, regions of both the 
incident and reflected fields. This theory has been applied 
to the study of two-dimensional aperture antennas by Clark for 
near zone scattering,^ Lysher used Babinet's principal to 

obtain the near zone fields in the rear hemisphere of a re- 

2 

f lector antenna by analyzing the fields of an aperture. In 
the aperture model, a plane wave, or in Lysher's case, a 
spherical wave, illuminates an opaque screen with an aperture. 
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The antenna beam is represented by the portion of the wave 
passinu: through the aperture. Clearly, this model only repre- 
sents one shadow boundary properly. Reflector antennas use a 
shaped reflector to focus the diverging rearward directed field 
of a feed into a collimated forward directed beam,'^ Two 
shadow boundaries are associated with this geometry, one the 
incident shadow boundary caused by the reflector shadowing the 
feed field, and the second, or reflection, shadow boundary 
between the reflected, focused, beam and the surrounding 
shadow. Note that the incident shadow boundary is improperly 
positioned in aperture' models of reflector antennas. Hatnasiri, 
Kouyoumjian and Pathak have applied the geometrical theory of 
diffraction to the study of the wide angle side lobes of re- 

4 

flector antennas in the far field. 



II* The Near Fields of Three-Dimensional ftef lector Antennas 
The analysis of a three-dimensional problem by the geo- 
metrical theory of diffraction is an extension of two-dimensional 
diffraction theory. The third dimension of the problem is 
introduced using either of two methods. First, if the field 
point is not near a caustic, the concepts of edge diffraction 
may be used. If the field point is near a caustic, the usual 
edge diffraction concept may be extended by use of an equivalent 
current technique. If the geometry involves a shadow boundary 
near the caustic, the equivalent current method yields the 
correct result only in the near half of the Fresnel zone. 

The general method used in this paper is to model the 
radiated field of the antenna by the geometrical optics field 
and the diffraction fields associated with the shadow boundaries 
of the geometrical optics field.. The diffraction fields of the 
antenna will be calculated from the position of the shadow 
boundaries, and these fields will be summed with the geometrical 

optics fields to obtain the total Fresnel zone fields, 

A, Basic Concepts 

I* Revie\y of Half-Plane Diffraction Mechanisms 
Figure r illustrates the various diffraction regions when 

t 

a point source illuminates a semi-infinite, planar, opaque, 
reflecting screen. For the point P in region I, the field at 
P is composed of a direct ray from S. and a ray reflected from 
the plane. In going from region I to region II, the reflection 
shadow boundary is crossed and now only the direct ray illumi- 
nates P. Finally going to region III, we cross the incident 
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Figure 1: Point Source Illuminating a Half Plane 

shadow boundary! and P is in the geometrical optics shadow of 
the half plane. In addition to the geometrical optics rays, 
there exist a class of rays caused by edge scattering, 

5 

Sommerfeld first solved the scattering or diffraction 
problem for a source distance of infinity or plane wave inci- 

g 

dence, Pathak and Kouyoumjian recently obtained an asymptotic 

« 

solution for a finite source distance, Tiie special case of 
their solution for illumination normal to the edge will be 
used in the present analysis, 

is the field of the incident scalar wave at 
the edge, the diffracted fields are given by 
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where the scalar diffraction coefficients are 
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and ^ f S, and S' are shown in Figure 1 F(ldLai;) is a 
correction factor that is needed only near the respective 
shadow boundary* When kLa^^lO, F(kLa'''} can be replaced by 
unity. 

The following coordinate transformations make the notation 
more compact and reinforce the concept that the diffraction 
fields are intimately related to the location of the shadow 
boundaries. The first term in (2) can be associated with the 
incident shadow boundary and the second with the reflection 
shadow boundary. We can relate these terms to the angular 
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distance from the shadow boundary. 



(G) 

4 + 4 * ’ 6 ^+*^ 

(7) 


Where 0^ , 0* arc the angles from the field point to the inci- 
dent and reflection shadov/ boundary, respectively, 0 is 
positive in the shadow region. The diffraction coefficient 
is separated into the incident and reflected components, which 
may be expressed as 


srBFiT 


where 


O*^ = 2. • 

From (2) and (8) we see 

D»CiV>o'-(60TD''(e“). 


( 8 ) 


(9) 


( 10 ) 


Kouyoumjian and Pathak have constructed a simple dyadic 
diffraction coefficient for vector fields using the soft and 
hard boundary scalar diffraction coefficients. 

where 

^ is the unit vector tangent to the edge. 
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A A ^ A 

p = c X I where I is the unit vector in 
the direction of incidence, and 

A A A 

pd = e X d where d is the unit vector in 
the direction of diffraction. 

Using the dyadic diffraction coefficient, one finds 






s(s+i‘) 


e ^ . 


( 12 ) 


To find the diffraction fields of a parabolic reflector, the 
above results must be extended to include the effect of a 
curved screen with a curved edge, 

2, Effect of Surface Curvature 
The diffraction fields act to eliminate the discontinuity 
at the shadow boundary. On the shadow boundary sufficiently 
removed from the edge, the diffracted field is exactly half 
the difference between the geometrical optics fields on the 
illuminated and shadow sides of the boundary. The diffracted 
fields at the shadow boundary of a half plane illuminated by 
a point source are given by (12) and the incident geometric 
fields on the shadow boundary by 




s' 


( 13 ) 


Since on the shadow boundary sufficiently removed from the edge 


then 



( 14 ) 
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(15) 


A. 

The last term in (15) is obtained by solvinf? (8) for a = 0. 

L is forced to depend on the source and observation distances 
to properly match the decay of the diffracted wave with dis- 
tance to the decay of the sceometrical optics wave. 

The key point is that L is determined by the divergence 
of the geometrical optics component. The divergence of the 
geometrical optics component at the incident shadow boundary 
is determined solely by the actual source and observation 
point distances. The divergence of the reflected field, 
however, depends on the curvature of the reflecting surface 
as well as the source and observation point distances. This 
effect of surface curvature can be included by determining 
the location of the image of the source and computing the di- 
vergence of the reflected ray using the image distance. Since 
the parameter L depends on the divergence of the associated 
geometrical optics field, L also depends on the image distance. 
Thus L will be different for the incident and reflection dif- 
fraction coefficients if the screen is curved. For the para- 
bolic dish with the feed at the focus, s* for the reflection 
diffraction coefficient will be infinity because the reflected 
geometrical optics field is that of' a plane wave, 

3» Effect of Edge Curvature 
One of the basic concepts in geometrical optics is the 

-9_ 


power conservation in a tube of rays. The assertion is since 
power flows only alon;; rays and not transverse to them, the 
power flux in a region bounded by rays is a constant. Thus 
the power and the fields at two points can be related to the 
cross section of the ray tube at these points. 



Figmre 2; Astigmatic Ray Tube 

Given the astigmatic ray tube shown in Figure 2 the ratios of 
the areas of the surfaces A, B may be expressed 

B 

Since the power flux is constant, the power density is inversely 
proportional to the area, and the field Intensity is inversely 
proportional to the square root of the area. 

Lines 1 and 2 in Figure 2 are referred to as caustics of the ray 
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tube, and t their associated caustic distances. 

In the diffraction nroblem the diffracting^ edrre is one 

caustic of the ray tube# The other caustic is determined by 

the i^eomctry of the edcc# A new formula for the second caus- 

7 

tic distance has been developed usinio: differential geometry# 

JL 

f ' ^ 




where 

j) is the caustic distance, 

1 is the distance from the diffracting point 
on the edge to the source, 

^ is the radius of curvature of the edge, 
n, e are the unit vectors normal and tangent 
to the edge, 

^ -t 

I, d are unit vectors in the direction of the 


incident and diffracted rays, and 



Applying (18) to a disc with source on the axis, one finds 
that the position of the caustic of the diffracted rays is the 
disc axis* 

Since the geometry is symmetrical about the z axis, one 
can choose any plane through the z axis* When the field point 



(19) 

( 20 ) 








where 


YfZ are the coordinates of the field point, 
A. is the radius of the disc, 
b is the nortiial distance from the disc to 
the source, and 

I is the angle between d and the z axis* 
From (18) the caustic distance p , '%hown in Fig. 3 is 

} ' r” 


L, -L' . 

0 R A 


Q » — ^ (24) 

j sii^ r • 

Thus is the distance on a line passing through the edge, P, 
and the z axis, from the edge to the z axis* 

From similar triangles, 

-4—1-IA- 

JiJ ' 1 y • 

In this equation J is negative if the caustic and P lie on the 
same side of the edge* Substituting this result into (l) with 
S* = y , we obtain 



The terra in the brackets provides for the ~ phase 


jump en- 


countered on passing throu$;h the caustic* It is included 
only if the caustic lies between the edge and P, 

8, E-Plane Analysis 

1, Heview of the Two-Point Technique 
0 

Keller developed an extension to Fermat's principle 
which relates the angle of diffraction to the angle of inci- 
dence. This extension states that a ray striking an edge ex- 
cites a fataily of diffracted rays which form a cone of half 
angle equal to the angle between the incident ray and the edge. 
Thus, for all field points not on a caustic of the edge, one 
can construct a finite number of edge diffracted rays which 
reach that field point. For the straight edge, only one ray 
fits this condition, while for the circular disc two rays 
are sufficient to describe the diffraction effects. 



(Rront View) 
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The fields will be calculated in the region y > 0, z > 0, 
Pattern symmetry allows these results to be reflected in the 
z axis to obtain the complete pattern. In the E-plane from 
Figure 4, the two points Q 2 lie on the y axis and e^, 

^ A 

are +x, -x respectively. The field point is at y, z; the 
subscript n identifies the diffracting point; and the super- 
scripts i, r associate the term with the incident or reflection 
shadow boundary. Since © is positive in the shadow region, 
and is contained in the yz plane, equation (25) reduces to 




•iX 


Using ( 10 ) to resolve into its components and (6) and 
( 7 ) to relate ^ *'nd O^we find 

- (-'T e;)] X 




where s ,© « *0 are obtained from the coordinates of the 

n n n 

field point using the following transformation and Figure 5 
0 a - AnclA/s/ ” 


'A 

The sign of© must be carefully preserved. 




- 15 - 




Fii^ure 5: Geometry of Two-Point E-Plane Analysis 
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^ is shown in 


Fii^ure 


2^ 

Kesolving P 


into axial and transverse components 


= ©AS 


S S\t^ . 


Thus 



(32) 


( 33 ) 


and 


e, (y,s)=-|. " D"( 

* AS»l 




( 34 ) 


2. Review of the Equivalent Current Technique 
When the field point is on or near the axis (corresponding 
to a caustic) the two-point method fails* For points in the 
vicinity of the axis, a method based on equivalent current 
flowing on the rim of the dish is available* 

The diffraction field of a half plane may be represented 
by equivalent currents flowing on the edge of the plane* To 
derive these currents consider the diffraction field when the 
electric field is parallel to the edge 





(35) 
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and when the magnetic field is parallel to the edge, 

-jAs 


( 36 ) 


These fields can be considered to be due to currents of elec- 
tricity and magnetism flowing along the edge* 




i 

) «}• 


^Wk. 


(37) 
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el 
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(38) 


In the general case one must take the component of the field 
tangent to the edge, then equating (35) with (37) and (3C) 


with (38) 



(39) 


(40) 


A standard radiation integral can now be used to find the 
fields due to these currents* 

In the present problem these currents flow on the edge of 
a circular disc so that the vector and scalar products in (39) 
and (40) simplify to trigonometric functions* The field of the 
feed is assumed to be constant around the edge of the antenna, 
i*e* the feed E-plane and H-plane patterns are identical* 

- The electric field due to the electric and magnetic currents 
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are given by the following integrals: 


-llt/V 




iii 

ifTf 




n (<!'■') A 4' A i>'. 


The above integrals assume that the field point is in the far 
field of the incremental sources even though it is in the near 
field of the totality of sources. As long as the incremental 
source to field point distance is greater than l.G wavelengths, 
this is a valid approximation. Substituting (39) and (40) 


in (41) and (42) gives 






-\U 


On the axis of the antenna the above integral simplifies 
considerably. The slant range, r, is a constant and the vector 


products reduce to 


A A 




j>* - 


A ct ^ a A > 
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Substituting; these results in the radiation integral and per- 
forming the integration, one obtains for the axial case 


pit, 

c - ^ 





( 46 ) 


The fields in the transition zone near the axis where the 
two-point method gives questionable results may be calculated 
using an extension of the technique used on the axis. We can 
compute the fields using the radiation integrals, inserting 
the actual incremental .source to field point distance in the 
integrals to account for the new location of the field point 



Figure 6: Determination of r Near z Axis 

Prom Figure 6 for a pattern in the yz plane 
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( 47 ) 



( 48 ) 
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( 49 ) 


For field points close to the axis D_(Q) is relatively constant, 

I 

We can take D_(Q) constant at its value for the on axis case 

§ 

Os (Q^ « (6‘) )> ( 50 ) 

H 

where e , O^are obtained from ( 30 , ( 31 )- for y = 0 , 

The equivalent current method will fail in the far field 
when a shadow boundary and a caustic of the diffracted rays 
are parallel. 

In the case of back scatter from the two-dimensional strip, 
one can show that the geometrical theory of diffraction pro- 
vides proper transition to the far field behavior; and the 
same is true of a three-dimensional problem with diverging 
shadow boundaries. 

For the strip, consider Figure 7. With the electric 
field parallel to the strip edge, the diffracted fields are 




( 51 ) 


since <1)*=^^ , s' = 00 for a plane wave at normal incidence. 
The reflected fields are 

£• V®' ^ ( 52 ) 
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A 

V 



Figiure 7; Conducting Strip Illuminated by a Plane Wave 
Referring to (2), (3) when s* and L = s 

z, 

f -it'. 1 

As we go to the far field on the axis 



(53) 



and the second 


“ *ri“ > 

term in predominates 


(54) 
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Fnrj small, the second term In (55) can be represented as 

'I Z «• 




When the diffracted fields from the two edges interact with 
the geometrical optics field, the first term in (58) just 
cancels the geometrical optics field and we are left with 

which is the expected result* 
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Now for the three-dimensional case with diverf^inp shadow 





Insertin.fl (61 ) into the radiation intecral 



♦'fre j cosTIii! ' 
r ^ 


Finally 



(fT^ 


( 62 ) 


(63) 


The result is correct to within a constant factor. The correct 
constant would be obtained if the magnetic current had also 
been included. 

When computing the field near the shadow boundary, the 
equivalent current concept fails since the diffracted fields 

f 

do not have the appropriate ~ dependence for the two- 

iTt 

dimensional geometry, and hence these diffracted fields cannot 
be equated to the fields of an infinite current source. 

HI, Near Fields of a Uniform Antenna 

In this section, the near fields of a uniformly illuminated 
antenna 10 wavelengths in diameter will be calculated. The 
near fields will be calculated for both the aperture model and 
the reflector geometry. Recall that the aperture model con- 
siders only the shadow boundary associated with the main beam, 
while the reflector geometry also includes the effect of the 
shadow of the reflector in the feed pattern.. Direct feed 
radiation and aperture blockage effects will not be considered. 
The equations developed in Section II were programmed in 
Fortran IV and run on the I®i 360-75 of the IHCC» The electric 
field was computed for planes 2, 10, 20 and 40 wavelengths in 
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front of the antenna usinr? the two-point method off the axis 
and the equivalent current method on the axis* The results 
are compared with similar data calculated by using plane 

wave spectrum concept and also aperture integration. Two 
results are presented for each distance, one including the 
effects of the incident shadow boundary and one neglecting it. 

The agreement between the diffraction theory results and 
the plane wave spectrum-aperture integration is very good. 
Resolution of the slight difference between the results would 
require resorting to more elegant methods of field analysis. 

The plane wave spectrum and aperture integration techniques 
both use the physical optics approximations to the actual 
field in the aperture. The assumption that the field is con- 
stant across the apertvire neglects the interaction of the 
aperture edges with the field. An exact analysis using an 
integral equation or eigenfunction technique would resolve 
the question. 

Figure '9 compares the GTD and aperture integration 
results at a range of 2 wavelengths from the antenna. The 
equivalent cxirrent technique was used to obtain field values 
for points out to one wavelength on either side of the axis. 

At this range, the two-point method and the equivalent current 
method agreed to within 0.2 A of the axis. In all the figures 
the erroneous two-point data near the axis is shown by the 
dotted lines. Figure !0 compares the GTD and aperture 
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integration results 10 wavelencths from the antenna, A.^ain 
the equivalent current tochnioue as^reed with the two-point 
method beyond 0,2 A from the axis. Note that the GTD solu- 
tion shows that the field on the axis is less at 10 X than at 
both 2 X and 20 X • Figures 11 and 12 show the GTD results 
for distances of 20 and 40 wavelengths respectively. At 20 X 
the two-point method began to show some error 0,5 A from the 
axis and got worse closer to the axial caustic. The error 
was appreciably greater than in the 2 X and 10 X cases. At 
40 \ the two-point method began to fail 1 X from the axial 
caustic. At this range the two-point method error near the 
axis was greater than in any of the smaller ranges. 

Figure 13 shows the field on the axis as a function of 
range. The parameter X is the ratio of the distance to the 
far field criterion. Note the oscillatory behavior of the 
field near the antenna. The field reaches the furthest 
maximum at 22 wavelengths and then should show ^ decay. The 
failure of the equivalent current method is clearly seen in 
this figure. Note the increasing field with distance beyond 
160 wavelengths. 

Agreement between the GTD and aperture integration re- 
sults is quite good. Reference to Tables I - IV will reveal 
that the effect of the incident shadow boundary is strongest 
at the greatest transverse distances ndiere the reflection 
shadow boundary term is small. 


- 27 - 



RELATIVE AMPLITUDE 



Figure 9: Electric Field on a Plane 2 Wavelengths 

From a Uniformly Illuminated Antenna 
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TABLE I 


NEAR FISLB OF UNIFORJ’f CIRCULAR ANTENNA ON A PLANE 


SURI'^ACE 2 '.VAVELENGTHS FROM THE ANTENNA* 


Transverse 

Distance 

Geometric Theory 
of Diffraction 

Plane Vtave 
Spectrum 

Aperture 

Intecrration 


RSB 

ISB + RSB 



0,0 


1.3420 

1,3101 

1.3356 

0.5 

.9290 

.8912 

.9186 

.9170 

1.0 

1,0470 

1.0710 

1.0547 

1.0640 

1.5 

.9732 

,9612 

,9679 

,9670 

2.0 

.9993 

.9986 

1.0022 

. 99o3 

2,5 

1 .0500 

1.0770 

1.0566 

1.0663 

3.0 

.8875 

,8387 

.8771 

,8666 

3.5 

1.1290 

1.1750 

1.1379 

1.1447 

4.0 

1.1020 

1.1430 

1.1066 

1.1083 

4.5 

.8073 

,8049 

.8102 

.8115 

5,0 

.4930 

.4420 

.4765 

.4751 

5,5 

.2921 

.2312 

.2833 

.2911 

6.0 

.1699 

,1030 

,1713 

,1756 

6.0 

.1020 

,0364 

,1060 

.1093 

7.0 

.0686 

.0115 

.0737 

.0771 

7.5 

.0511 

.0061 

.0569 

.0602 

8.0 

,0403 

,0034 

,0457 

.0493 

8.5 

.0326 

,0044 

.0373 

.0411 

9.0 

.0269 

.0055 

.0313 

.0346 

9.5 

.0245 

.0065 

.0260 

.0293 


* Data and curves show GTD equivalent current results for 
the first wavelength from the axis, and two point method 
results for the remainder of the plane* 
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RELATIVE AMPLITUDE (''O^'^VmETER) 



■ 6 -4 -2 0 Z t 

TRAflSVERSE DISTANCE ( WAVELENGTHS) 


Figure 10; Sleetric Field on a Plane 10 Wavelengths 


From a Uniformly Illuminated Antenna 



TABLE II 


NEAR FIELD OF UNIFORM CIRCULAR ANTENNA ON A PLANS 
SURFACE 10 WAVELENGTHS FROM THE ANTENNA 


Transverse Geometric Theory 
Distance of Diffraction 


Plane Wave Aperture 
Spectrum Inte<rration 


RSB ISB + RSB 


0,0 


.9560 

1.0212 

1.0302 

0,5 

.8642 

.8827 

.8843 

.8873 

1.0 

.9963 

1.0090 

,9815 

.9837 

1.5 

1.0910 

1.1180 

1.0879 

1.0902 

2.0 

1.1500 

1.1620 

1.1468 

1.1470 

2.5 

1.2830 

1.3430 

1.2978 

1.3032 

3.0 

1.1070 

1.1350 

1.1293 

1.1308 

3.5 

.8655 

,8556 

,8817 

,8780 

4.0 

.7853 

.7963 

.8008 

.8011 

4.5 

.6207 

.6093 

,6138 

.6143 

5.0 

.4611 

.4242 

.4244 

.4254 

5.5 

.3995 

.3843 

.3836 

,3860 

6.0 

.3151 

.2964 

,3050 

.3063 

6,5 

.2044 

.1658 

.1954 

.1931 

7.0 

.1785 

,1527 

.1770 

.1773 

7.5 

.1714 

.1549 

.1733 

.1763 

8.0 

.1261 

,1023 

.1259 

.1283 

8.5 

.0829 

,0481 

.0811 

.0812 

9.0 

.0830 

.0601 

,0845 

.0850 

9.5 

.0893 

,0726 

,0927 

,0946 
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TRANSVERSE DISTANCE ( WAVELENGTHS ) 


Figure 11: Electric Field on a Plane 20 Wavelengths 

From a Uniformly Illuminated Antenna 



TABLE III 


NEAR FIELD OP UNIFORM CIRCULAR ANTENNA ON A PLANE 
SURFACE 20 WAVELENGTHS PROM THE APERTURE 


Transverse 

Distance 

Geometric Theory 
of Diffraction 

. Plane V/ave 
Spectrum 

Aperture 

Integration 

0.0 

RSB 

ISB + RSB 
1.8500 

1.8419 

1.8472 

0.5 

1.7650 

1.8770 

1.7179 

1.7226 

1.0 

1.4030 

1.4550 

1.4220 

1.4242 

1.5 

1.1260 

1.1310 

1.1528 

1.1527 

2.0 

1.0440 

1.0510 

1.0769 

1.0768 

2.5 

1,0560 

1.0840 

1.0960 

1.0975 

o.O 

.9909 

1.0200 

1.0282 

1.0305 

3.5 

.8119 

.8197 

.8355 

.8570 

4.0 

,6049 

.5843 

.6069 

.6064 

4.5 

,4930 

.4641 

.4757 

,4741 

5.0 

,4983 

,4866 

.4586 

.4581 

5.5 

.4618 

.4600 

.4346 

.4352 

6.0 

.3716 

,3633 

,3497 

.3500 

6.5 

.2332 

.2312 

. ,2380 

,2372 

7.0 

.1862 

.1559 

.1773 

.1765 

7.5 

.1916 

.1768 

.1897 

.1907 

8.0 

.1989 

.1922 

.1983 

,2004 

8.5 

.1710 

.1624 

.1694 

.1712 

9.0 

.1189 

.1023 

.1160 

.1167 

9.5 

.0781 

.0512 

,0758 

.0749 
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Figure 12: Electric Field on a Plane 40 Wavelengths 

FVom a Uniformly Illuminated Antenna 
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TABLE IV 


NEAR FIELD OF A UNIFORM CIRCULAR ANTENNA ON A PLANE 
SURFACE 40 WAVELENGTHS FROM THE APERTURE 


Transverse 

Distance 

Geometric Theory 
of Diffraction 

Plane Wave 
Spectrum 

Aperture 

Intee:ration 


RSB 

ISB + RSB 



0.0 


1.5000 

1.6528 

1,6536 

0,5 

1,6520 

1.4570 

1.6209 

1,6224 

1.0 

1.4630 

1.3940 

1.5304 

1.5315 

1.5 

1,3040 

1.3480 

1.3878 

1.3891 

2.0 

1.1320 

1.1620 

1.2071 

1.2080 

2.5 

,9502 

,9651 

1.0040 

1.0048 

3.0 

.7714 

,7718 

.7987 

.7992 

3.5 

.6121 

.5994 

.6131 

.6135 

4.0 

.4913 

.4689 

.4721 

.4721 

4.5 

.4274 

.4026 

.3947 

.3946 

5.0 

.4234 

.4046 

.3754 

.3754 

5,5 

.4085 

.3991 

,3811 

.3814 

6.0 

.4012 

• 39oo 

.3815 

.3820 

6.5 

.3824 

.3820 

,3625 

.3631 

7.0 

.3404 

.3391 

.3222 

.3228 

7.5 

,2825 

,2779 

,2658 

.2661 

8.0 

.2169 

.2071 

.2020 

.2020 

8,5 

,1550 

.1388 

.1433 

.1427 

9.0 

.1132 

.0918 

,1067 

.1058 

9.5 

.1058 

,0886 

,1054 

.1046 
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RELATIVE AMPLITUDE (<I8) 



AXIAL DISTANCE- 2 

Figure 13 : Electric Fields on the Axis of a 

Uniformly Illuminated Antenna 
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IV, Near Fields of an Actual Antenna 

The near zone fields of an actual parabolic antenna were 
computed and compared with measured results. The antenna 
studied was a 24-inch diameter 8-inch focal lencth spun alu- 
minum dish with sharp edges. The antenna was fed with a 
section of open ended wave guide with a plain flange. The 
feed was supported on a tripod constructed of 3/8-inch dia- 
meter polystyrene rods. The rods were fastened to the dish 
about 2 inches from the edge to reduce interaction with dif- 
fraction effects. Feed patterns are presented in Figure 14 : 

The antenna was driven by a Sylvania Model 3200 avalanche 
diode oscillator mounted on the feed. 

Direct exurrent was fed to the diode oscillator through a 
miniature coaxial cable. The avalanche diode power source 
was chosen to minimize the cross section of components in the 
main beam of the antenna with the intent of minimizing spurious 
pattern disturbances due to aperture blockage. The antenna 
measurements were made in the Student Antenna Laboratory, 

Room 731 in the Electronics Laboratory, This room occupies 
two floors of the building providing a 20 x 20 x 40 foot 
volume for antenna measurements. The room surfaces are con- 
ventionally constructed with asphalt tile floor, plaster walls 
and a suspended grid ceiling, Eccosorb has been applied to a 
six-foot square area of the side walls and one end wall to 
reduce specular reflections. Facilities available include a 
pattern recorder for plotting polar antenna patterns. This 
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RELATIVE AMPLITUDE (dB) 



Figure 14; Open End Waveguide Peed Pattern 
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plotter was adapted to allow plottine straight line scans. 

The antenna was fHOunted on a lO-foot mast in the center of 
the room. The probe consisted of an open-ended wave iruide 
section mounted on a moveable carriacre and coupled to a cry- 
stal detector. Adequate signal level was obtained at all 
ranges. 

The development of Section II must be modified to in- 
clude the illumination taper introduced by the feed pattern. 
Two corrections are required, one is that the field at the 
reflector rim is reduced, and the second is that the geomet- 
rical optics field must vary across the aperture. The re- 
duced value of the field at the edee can be introduced by a 
constant multiplicative factor in the calculations. 

To find the geometrical optics field, we need to find 
the power density in the aperture as a function of distance 
from the axis, Tlie power density gives the field in the 
aperture, and we can then use the principal of stationary 
phase to assert that the seometrical optics field can be 
related to the field in the aperture at the stationary point. 
The stationary phase arguement that follows may be used to 
derive (17), 

The radiation integral is of the form 



- f\ 
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For the stationary phase principle to apply 


3 ^ 


-0 , 


( 65 ) 


must he satisfied at some point called the stationary 
point. In our problem 




( 66 ) 


where 



in order that stationary phase concept be applicable. 


(67) 


V « j * • (68) 

This phase function fulfills (65) for y^ = y. The essence of 
the stationary phase arp:«ement is that the limits of (64) can 
be reduced so that only an elemental area about is considered# 
This is valid because the region where (|) is slowly varying 
contributes more to the integral than do the regions where ^ 
is varying rapidly. Once the range of integration is restricted, 
F(y*,0) can be considered constant and taken outside the inte- 
gral; and (64) becomes 


M- ('5 .*■) = ’•1'^) Aj'. 


(69) 
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Due to the stationary process the value of the integral will 
not change significantly when we let the limits of integration 
go to infinity. 

oo ^ ^ 

^ 

J (70) 

-««a 

The definite integral is easily evaluated. Thus we see that 
the field depends on the forcing function at the stationary 
point. 

To find the aperture field, one must first find the 
angle ^ in Figure 15 use the far field pattern of the feed 
to find a reference field at each value of the angle; and then 
account for the ray tube spreading over the various distances 
from the feed to the reflector. Since the parabola collimates 
the diverging bundle of rays from the feed, there is no field 
decay in the portion of the path from the reflector to the 
aperture. Thus the povrer density in the aperture is the same 
as that at the reflector. The angle iff can be simply related 
to the distance of the parallel ray from the axis and the 
focal length of the reflector. 






zi 


( 71 ) 


where y is one coordinate of the field point, and 
^ is the reflector focal length. 

Since we know the far field pattern of the feed, one can use 
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curve fitting techniques to obtain an analytic function for 
the portion of the pattern illuminating the reflector. In 
the E-plane, this function is 

( 72 ) 

Since is unity at ^ = O, it is convenient to normalize the 

field in terms of the field in the center of the aperture. 

Thus represents the magnitude of the field on a spher- 

ical surface whose radius is the focal length and centered 
on the feed. This surface is represented by the dotted circle 
in Figure 15. The decay in power density as the ray travels 
from the surface to the reflector may be found using power 
conservation 





(73) 


where is the amplitude at the reflector, 

LKg is the amplitude on the surface, 

5^/v are shown in Figure 15. 

/i. may be found from simple geometry to be related to the dis- 

« 

placement of the reflected ray. 





(74) 


These expressions may be combined to give the field in 
the aperture as a function of the distance of the aperture 
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point from the z nxis 




■k 


Z rti?CT/»w 



( 75 ) 


This expression for the geometrical optics component was summed 
with ( 35 ) and (48) to give the total near fields at 16,6, 37,6, 
74, 117,5 and 188 wavelengths. These distances were chosen 
from Figure 16, a plot of the measured fields on the axis 
versus distance from the antenna. The distances represent 
maxima and minima in power density. Comparisons are shown 
in Figures 17 through 21, 

The agreement between measured and computed results is not 
as good as one would have expected. Direct feed radiation and 
feed aperture blockage were first suspected of causing the 
discrepancy. The magnitude of feed radiation was estimated 
using geometrical optics techniques. Since the worst case 
radiation was of interest, it was assumed that the feed pattern 
was uniform at the rear axial level over the interesting por- 
tion of the rear hemisphere. Since the feed field at the 
center of the reflector was normalized to one, the field at 
any point is 

where ^ is the distance from the feed to the 
field point in wavelengths. 
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Figure 16 , Measured Axial Power Density for an Antenna 

With Tapered Illumination ^ measured data in Following 
Figures has been Normalized to get Best Fit in the Near 

Zone. 



5 


0 


-5 



TRANSVERSE DISTANCE— (WAVELENGTHS) 

Figure 17: Electric Field on a Plane 16«6 Wavelengths 

Froni an Antenna with Tapered Illumination 
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RELATIVE AMPLITUDE (dB) 



TRANSVERSE DISTANCE (WAVELENGTHS) 

Figure 18; Electric Field on a Plane 37.6 Wavelengths 
From an Antenna with Tapered Illumination 
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Pigvire 19: Electric Field on a Plane 74 Wavelengths 

From an Antenna with Tapered Illumination 
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RELATIVE AMPLITUDE (dB) 



Figure 20: Electric Field on a Plane 117,5 Wavelengths 

From an Anteima with Tapered Illumination 
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RELATIVE AMPLITUDE ( dB ) 



TRANSVERSE DISTANCE ( WAVELENGTHS ) 


Fi^^ure 21: Electric Field on a Plane 188.0 Wavelengths 
Prom an Antenna with Tapered Illumination 
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For aspects close to the axis 


%- A , 


( 77 ) 


where d is the distance from the feed to 
the aperture. 

The results of the estimation are shown in Table V, 

To estimate aperture blockacre, the blocking components 
were conceptually replaced with an equivalent source. The 
source is planar, equal in size to the projected size of the 
blocking components and driven by fields equal in magnitude 
but opposite in phase to the fields present in the plane of 
the blocking components in the absence of the blocking com- 
ponents. The blocking components are the metal waveguide 
components and the polystyrene mount. This complex assembly 
was modelled by a 2 A diameter disc. Since the disc is small, 
the far field approximation is used to find its field. Again 
being interested in the worst case, the scattering pattern was 
assumed constant at its maximutn level. Use the same notation 
as for. the feed radiation 



Comparing the results of the above estimates with the 
diffracted fields in Table V, we see that aperture blocicage 
dominates over feed radiation near the axis. At the smaller 
ranges, aperture blocking is sufficiently strong to cause 
errors. However the discrepancy at 188 X shown in Figure 21 
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TABLE V 


COMPARISON OF FEED, APERTURE BLOCKAGE 
AND DIFFRACTION FIELDS 


Distance 

Feed 

Aperture 

Diffraction 

sA 

Wavelencrths 

Radiation 

B1 ockacre 

On Axis 

Outside ! 

16.6 

.045 

,270 

.2929 

,0384 

23,3 

.031 

.170 

,3525 

.0396 

37,6 

.018 

.100 

.4136 

,0584 

72.0 

.009 

.051 

.4892 

.0839 

117,5 

,005 

.028 

,5563 

.0589 

118.0 

.0036 

.018 

.6311 

.0925 



cannot be explained by this mechanism. The side lobe error is 
especially severe. 

At this point it was decided to calculate the expected 
pattern by a second method. The pattern at 188 X was chosen 
for comparison because it is well into the Fresnel zone, and 
Loinmel functions may be used to perform the radiation inte- 
gration in closed form. The results of a hand calculation of 
two points on the side lobe are sho^m in Figure 21. This cal- 
culation was based on a uniformly illuminated aperture. The 
results were modified to correspond with the GTD model for the 
tapered aperture by reducing the unifoimi illumination in the 
Lomroel calculations to the same level as the edge illumination 
in the GTD model. Since the GTD model does not include the 
effects of the illumination taper at the edge, this is a valid 
step. The Lommel results agree precisely with the GTD results 
showing that the experimental data are faulty. 

The experimental method was reviewed to determine the 
source of the deviation* After considering room reflections 
and plotter enrors, it became apparent through repeated study 
of the patterns that the 188 X pattern was characteristic of 
a far field pattern. The null between the main lobe and the 
first side lobe was especially strange since the distance in- 
volved was about one quarter of the far field criterion. The 
appearance of a far field pattern in the near zone brought to 
mind the antenna focusing procedure* Since space was restricted 
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in the laboratory and the location of the test antenna was 
fixed to stabilize reflection, the antenna was focused with 
the test antenna and the probe spaced by 18SA • The feed 
was positioned for best gain and side lobe ratio. This results 
in the far field pattern appearing in the Fresnel zone. It is 
interesting to calculate the effect of this feed position. 

Assuming that the feed position causes a phase taper in 
the aperture which is the best fit to a spherical wave with 
radius equal to the spacing between the test antenna and the 
probe, the maximum phase change across the aperture is given 
by 

> (79) 

where z is the distance between the test 

antenna and the probe, and 
A 

^ s 2 arc tan ’** , the angle subtended 
z 

by the test antenna at the probe. 

For a 19.1 X antenna focused at 188 \ , the phase taper is 
88^ The feed position for this phase error may be calculated 
by matching the phase shift for a central ray and an edge ray 
with the required error. To match an 88° phase shift, the 
feed must be displaced from the focal point of the reflector 
by .3 ^ or very nearly one centimeter. To check this result, 
use the thin lens equation 
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(80) 



where i is the ima.'ice (probe) distance, and 
O is the object (feed) distance. 

Using this relation, the required displacement of the feed 
is 0.22 A which agrees well with the results of phase taper 
matching. 

Thus, this antenna adjusted for best pattern at 188 
exhibits sufficient phase taper to give strong disagreement 
between measured results and computed results which assume 
no phase taper. This phase taper affects the patterns at 
all distances and is considered to be the mechanism that 
causes the disagreement bettreen the computed and measured 
data. 

V., Conclusions 

The geometrical theory of diffraction is an effective 
method of calculating the near fields of reflector antennas. 
This work extends the results of others to treat the forward 
axial region by including the reflected field of the antenna 
beam. The methods developed in this paper fail as the field 
point moves into the far zone. This is not a major problem 
because other geometrical theory of diffraction techniques 
are available there. 

Comparison of the results for a uniformly illuminated 
antenna with data obtained by aperture integration and the 
plane wave spectrum technique show that equivalent results are 
obtainable with an expenditure of much less computer time. 
These results prove that the geometrical theory of diffraction 
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is applicable in the near zone of reflector antennas. 

Since the phase characteristic of the antenna used to 
gather experimental data did not match the model used in the 
computation, the results for an antenna with tapered illumi- 
nation cannot be completely verified. The technique for com- 
puting fields for a tapered aperture should be improved using 
the traveling wave concepts used in the GTD analysis of the 
H-plane pattern of an open-ended rectangular wave guide. The 
experiment should be repeated using an antenna focused in the 
far field. 
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